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ABSTRACT 

Fast-spinning strongly magnetized newborn neutron stars (NSs), including nascent magnetars, are 
popularly implemented as the engine of luminous stellar explosions. Here, we consider the scenario 
that they power various stripped-envelope (SE) supernovae (SNe), not only super-luminous SNe Ic 
but also broad-line (BL) SNe Ibc and possibly some ordinary supernovae Ibc. This scenario is also 
motivated by the hypothesis that Galactic magnetars largely originate from fast-spinning NSs as 
remnants of SE SNe. By consistently modeling the energy injection from magnetized wind and ®®Ni 
decay, we show that proto-NSs with > 10 ms rotation and poloidal magnetic field of i?dip ^ 5 x 10^^ G 
can be harbored in ordinary SNe Ibc. On the other hand, millisecond proto-NSs can solely power BL 
SNe Ibc if they are born with Hdip ^ 5 x 10^"^ G, and superluminous SNe Ic with Hdip ^ 10^^ G. Then, 
we study how multi-messenger emission can be used to discriminate such pulsar-driven SN models 
from other competitive scenarios. First, high-energy X-ray and gamma-ray emission from embryonic 
pulsar wind nebulae can probe the underlying newborn pulsar. Follow-up observations of SE SNe 
using NuSTAR ^ 50— 100 days after the explosion is strongly encouraged for nearby objects. We also 
discuss possible effects of gravitational waves (GWs) on the spin-down of proto-NSs. If millisecond 
proto-NSs with Bdip < a few x 10^^ G emit GWs through, e.g., non-axisymmetric rotation deformed 
by the inner toroidal fields of Bt > 10^® G, the GW signal can be detectable from ordinary SNe Ibc 
in the Virgo cluster by Advanced LIGO, Advanced Virgo, and KAGRA. 


1. INTRODUCTION 

Time-domain astronomy is rapidly expanding. Var¬ 
ious transients are now being efficiently detected and 
newly discover ed by survey facilities . Swift (iGehrels et al.l 
l2004li , Fermi (lAtwood et al.l I20 09H . the Palomar Tran¬ 
sient Factory (PTF: iLaw et ijiLll2009tl and the Panoramic 
Survey T elescope fc Rapid Response System (Pan- 
STARRS: iHodapp et al.l I2004D . Deeper follow-up ob¬ 
servations are well functioning i n rad io to gamma- 
ray bands (e.g., iGreiner et al.l 120081 : iPerlev et al.l 
120111 : lAbevsekara et al.l 120121 1. Mo reover, Che r enkcw 
detectors like Super-Kamiokande (llkeda et al.l 1200711 
and IceCube (jAartsen et al.l 1201,511 are now gaz¬ 
ing at the neutrino sky, and the installation of 
gravitat ional-wave (GW) interferometers like Advance d 
LIGO ((Harry k, LIGO S ci entific Collaborationj 1201011 . 
Advan ced Virgo (lAccadia et al.l I2011I L md KA¬ 
GRA (|Somivall2012t l will be completed soon. These new 
and upgraded facilities can potentially probe hidden en¬ 
gines of violent astrophysical phenomena. 

Some targets among the promising sources of GWs 
and neutrinos are fast-spinning strongly magne- 
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tized proto-neutron stars (NSs) formed in collapsing 
stars. They have been proposed as the engines of 
luminous tran s ients, e.g., long gamma-ray bursts (L- 


1998 

; iZhang & MeszarosI 

20011: IThomnson et al. 

2004 

; IMetzger et al.l 120071: 

Bucciantini et al. 1200911 


includin g their sub-c la ss w ith low luminosity (LL- 
GRBs: [ Mazzali et al.l 120061 : iSoderberg et al.l 120061 
Toma et al.l 12(107 ), br o ad-line type Ibc SNe (BL-SN Ibc: 
Wheeler et al.l l2000l iThomnson et al.l 12004 iWooslevI 
201C), hydrogen-poor super l uminous SNe (SL-SNe 


Ic: 


--- .Kasen fc: Bildsted 120101 iPastorello et al.l 120101 

Quiinbv et al.l 12011 : liiserra et al.l l2013l~Nicholl et al l 

201,811 ('see also lMetzger et al.l 1201,^ iWang et al.l 120151) . 
The basic picture is that the rotational energy of 
proto-NS is extracted by the unipolar induction as 
magnetized wind or jet and is later dissipated by a 
physical process that has still to be constrained, result¬ 
ing in luminous electromagnetic radiation. However, 
no observational finding has been able to conclusively 
validate the pulsar-driven scenario so far. The question 
is how to discriminate newborn pulsar engines for 
each type of transient by using ongoing and upcoming 
multi-messenger observations. 

In this paper, using a semi-analytical model shown 
in the Appendix, we consistently calculate the multi¬ 
messenger counterparts from fast-spinning strongly mag¬ 
netized proto-NSs, focusing on the cases accompanied by 
stripped-envelope (SE) SNe. In the next section, we dis¬ 
cuss another important motivation of our study: the pos¬ 
sible connection between Galactic magnetars and pulsar- 
driven SE SNe. Then, we consider the SN counterpart 
and derive the parameter range of the pulsar-driven SN 
model consistent with the observed SN Ibc, BL-SN Ibc 
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Fig. 1. — All-sky event rate of explosive phenomena poten¬ 
tially explained by pulsar-driven SN mo del. The core- c ollaps e 
supernova (CC- SN) rate is obtaine d from [HorjuchLeL-MJ (12^^). 
We also refer to [Smit het__^ Il2011l'l. IG uetta &; Della Valid lj2QQ7i 'l. 
IWanderman &; PiranI 11201011 and ibiuimbv et al.l Il20l3l for the rates 
of supernova Ibc and the peculiar sub-class with broad-line features 
(BL-SN Ibc), low-luminosity GRB (LL-GRB), long GRB (L-GRB), 
and hydrogen-poor superluminous supernova (SL-SN Ic), respec¬ 
tively. The vertical dashed line shows the mean distance to the 
Virgo cluster (= 16.5 Mpc). 

and SLSN-Ic (Sec. 3)0 In Sec. 4, we show the de¬ 
tectability of the multi-messenger counterparts including 
the pulsar wind nebular (PWN) emission, GW emission, 
and neutrino emission. Based on the results, we discuss 
observational strategies for the multi-messenger search of 
fast-spinning newborn NSs in SE SNe and possible sci¬ 
entific impacts in Sec. 5. We summarize our paper in 
Sec. 6. 

2. CONNECTION BETWEEN GALACTIC MAGNETARS 
AND STRIPPED-ENVELOPE SUPERNOVAE? 

Confirming pulsar-driven scenarios is important in 
terms of understanding the origin of Galactic magne- 
tars. In the classical picture ({Duncan fc Thompson|[l9^ 
iThompson fc Duncanlll99§l . the magnetic field amplifi¬ 
cation is attributed to the proto-NS convection coupled 
with a differential rotation less than a few ms. Even in 
the absence of such rapid rotation, the magnetic fields 
could be amplified by the magnetorotational instabil¬ 
ity (e.g.. iBalbus fc HawlevI 1199^ lAkivama et al.l 120031 
iThompson et al.l 120051 iMosta et al.l 1201511 . The total 
magnetic field energy of a magnetar is estimated to be: 

Eb - BtySn X 47ri?„eV3 - 2.9 x 10^® erg , 

( 1 ) 

while the free rotational energy stored in the proto-NS 
is: 

frot.i ~ /(2^/Pi)V2 ^ 3.1 X 10^® ( erg. (2) 

Here, Bt is the inner toroidal field strength, I ^ 
1.4 X 10^^ g cm^ is th e momentum of iner¬ 
tia (jLattimer fc Prakashll20Qlh PI and Pi is the initial spin 

^ In this paper, we do not include SNe Ilb, which are usually 
listed as part of the SE SN family. 

® We take a fiducial NS mass and radius of Mns = 1.4 Mq and 
Rns = 12 km, respectively. 


period. Even cases with P{ > 10 ms have a free energy 
of 

^rot.i 105O 

erg, which is sufficient to power SN ex¬ 
plosions. 

Formation of fast-spinning strongly magnetized proto- 
NSs may not be rare. Population synthesis calcu¬ 
lations of Galactic NS pulsars showed that the ini¬ 
tial spin distribution is a Gaussian with a peak 
at ^ 2 00-300 ms and standard dev i ation of ^ 

100 ms (iFaucher-Giguere fc Kaspil 120061 iPopov et al.l 
l2fnnli . If this applies even outside the Galaxy then the 
formation rate of proto-NSs with Pi < 30 ms is: 

B-Pi <30 ms ^ 0-1 X "^CC-SN- (3) 


Here, 7 ?.cc-SN 0-03lo;o26 ^ is the core¬ 

collapse SN rate le.g.. lAdams et al1l2013ll . On the other 
hand, the formation rate of Galactic magnetar has been 
esti mated from the observed spin-down rate and remnant 
age (iKeane fc Krameill2008ll : 


"^magnetar ^ 0.1 X PcC — gN. (4) 


Although uncertainties in both Eq. ([3]) and (|1]) are fairly 
large, the formation rate of fast-spinning proto-NSs and 
Galactic magnetars is still consistent with the dynamo 
scenario. At this stage, there is no obser vational support 
on th is interesting possibility (see, e.g., iVink fc Kuiped 

[2onl . 


The progenitors of Galactic magnetars are consid¬ 
ered to be very massive stars with Mzams ^ 30 — 
4 OM 0 based on the fact that the y are obser v ed in 
young massive star clusters (e.g ., iFiger et al.l 120051: 
Gaensler et al.l 120051 iMuno et al.l 120061 iBibbv et ahl 


20081 iDavies et all 200S), and dis tributed in low Galactic 


latitudes ({Olausen fc Kaspill2014f) . Note that the fraction 
of massive stars with Mzams ^ AOMq is roughly ^ 10% 
of that of Mzams ^ SMq given the Salpeter initial-mass 
function. Such massive progenitors with approximately 
a solar metallicity are considered to evolve into Wolf- 
Rayet stars (WRs ), and end their lives as SE SNe (e.g., 
iHeger et al.l 1200311 0 The observed event rate of ener¬ 
getic SE-SNe and associated high-energy transients are 
relatively low: 


■^^BL-SN Ibc 0.04 X Pcc-SN, (5) 

■T^LL-GRB 0.01 X Pcc-SN, (6) 

■T^L-GRB ~ 'T^SL-SN Ic lO”'^ X PcC-SN, (7) 

(IGuetta fc Della Valid 120071 IWanderman fc Piranl[201(11 

iSmith et al.ll20111lQuimbv et al.ll2013[l . As for L-GRB, a 
jet be aming factor of ^ 100 is assumed (e.g. IGuetta et al.l 
1200511 . Even if these transients are powered by newborn 
magnetars, they explain only a minor fraction of the total 
magnetar abundance. On the other hand, ordinary SN 
Ibc can meet the magnetar formation rate as high as Eq. 
( 0 : 

72-ibc 0.26 X TT-cc-sn- (8) 

Actually, e.g., iMaeda et al.l (|2007ll proposed a newborn 
magnetar as a relevant energy source of a type Ibc SN 
2005bf. In this regard, it is important to show in what 


^ We note that a significant fraction of the observed SE 
SNe may be from close binary systems fe.g.. [Eldridge et al.ll20^ : 
ISmith et al.ll20Tiri . 
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Fig. 2. — SN light curves of pulsar-driven SN models with ms 
rotation and different magnetic field strength. The upper and lower 
gray lines show the observed SL-SN-Ic PTF 09cnd and BL-SNe Ic 
1998bw, respectively. 


parameter range ordinary SN Ibc is compatible with the 
pulsar-driven scenario and how to identify the underlying 
newborn pulsars. 

We should note that it has been pointed out that 
pulsar-driven models cannot reproduce the observed light 
curves of SNe Ibc and BL-SNe Ibc, in particular the 
late-time behavior > 100 days after the explosion (e.g., 
iSollerman et al.ll200^HiTserra et aI1l2013D . In such a late 
phase, however, theoretical modeling of the optical light 
curve is still largely uncertain (see Sec. I3.2l and Sec. 15.21) . 
In this paper, we focus on the optical light curves around 
the peak where the diffusion approximation is robust. On 
the other hand, in the late phase, the SN ejecta becomes 
almost transparent for nascent PWN emissions in hard 
X rays and gamma rays. They can be good probes of the 
properties of the underlying pulsar and one should take 
into account the fact that only a fraction of the energy 
can be converted into optical emission since high-energy 
emission escapes. 


3. PULSAR-DRIVEN SUPERNOVA SCENARIOS 
As popularly discussed in the literature (e.g. 


Ostriker & GunnI 119691 

Thomoson et al.l 12004 IWooslevI 

20101: iKasen & Bildsten 

120101 IWansi et al.l 120151). verv 


bright SNe could be explained by the pulsar-driven SN 
model with P; less than a few ms. In this scenario, the 
peak luminosity of the pulsar-d riven SN can be estimate d 
as LPf « (|Kasen fc Bildsteiill2010D . 

or 

xf__") f 1^9) 

\2 X 10® cm s~^ J \0.2 g“^ cm^y 

Here, 


j.em 

^sd 





( 10 ) 


is the dipole spin-down timescalej^ and 


Kic ^ 20 days 


jy±_ 

5 Mq 


1/2 






2 X 10® cm s“^ 
1/2 


- 1/2 


0.2 g“i cm- 


( 11 ) 


is the photon diffusion time from the ejecta. 

In this work, we numerically calculate light curves of 
SNe driven by fast-spinning strongly magnetized new¬ 
born NSs embedded in SE progenitors. Details of the 
model description are given in the Appendix. We as¬ 
sume that the energy injection is caused by spherical 
winds rather than jets and both ®®Ni decay and mag¬ 
netized wind are taken into account as energy sources. 
The thermalization of the non-thermal emission is ap¬ 
proximately taken into account, and the optical SN emis¬ 
sion and early non-thermal nebular emission are obtained 
consistently. The effect of GW spin-down is incorpo¬ 
rated in a simple parametric fo rm. The present model 
is based on iMurase et al.l (|2015f) but with several refine¬ 
ments, e.g., including the effect of ®®Ni decay. The simple 
model allows us to explore a wide parameter range; the 
initial spin of Pi = 1 — 30 ms, poloidal magnetic field of 
Bdip = G, SN ejecta mass of Mej = 1 — 10 Mq, 

^®Ni mass of = 0.05 — 1.0 Mq, SN explosion 

energy of Sgn = 10®^“®® erg, and graybody opacity 
Kt = 0.05 - 0.2 g-i cm®. Note that Kj 0.1 and 
0.2 g“^ cm® corresponds to electron scattering for singly 
ionized and fully ionized helium, respectively, and can be 
smaller for, e.g., a partially ionized C- or 0-dominated 
ejecta. 

In Fig. [21 we show some light curve exam¬ 
ples of the millisecond-pulsar-driven SN model. The 
thicker red lines correspond to larger magnetic fields. 
The g ray lines indicate the observed SL -SN-Ic PTF 
09cnd (iQuimbv et al.ll2011l: lGal-Yamll2012l) and BL-SNe 
Ic 1998bw (jGalama et al.lll998ll . Pulsar-driven SNe be¬ 
come as bright as SL-SNe with Pdip ^ 10^"^ G. In such 
cases, a significant fraction of the spin-down luminosity 
needs to be converted into SN radiation. The pulsar- 
driven SN model with Pdip ^ 10^^ G and Pi 1 ms can 
reproduce the observed light curve of this event. 

We also consider the stronger case, where energy in¬ 
jection from fast-spinning NSs contributes to some BL- 
SNe Ibc and possibly ordinary SNe Ibc. For a hxed ini¬ 
tial spin, the peak luminosity becomes smaller with a 
stronger magnetic field since f®™ becomes smaller (see 
Fig. 121 and Eqs. I MTU)) . The physical reason is that 
the proto-NS spins down long before the photon diffu¬ 
sion time, and the injected energy by the pulsar wind 
is lost via adiabatic cooling. This means, on the other 
hand, that the injected energy is used for acceleration 
of the ejecta rather than SN radiation. Interestingly, as 
for Pi about a few ms and Pdip ^ 5 x 10^“^ G, the peak 
luminosity becomes pP®’’ < erg s“^ and the mean 


As for the spin-down luminosity, we use a formula motivated 
by up-to-date MHD simulations, which give a factor 9/2 l arge r 
value on average than the classical dipole formula (see Eq. IA2II . 
As a result, and Lsn become smaller by the same factor for 
a given and P[. This difference may affect the estimation of 

these parameters from observations. 
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Days since Maximum 


Fig. 3.— Comparison of observed SN Ibc and BL-SN Ibc (gray 
lines) and theoretical light curves. The red solid and dotted-dashed 
lines show wind-powered cases, whereas the blue dashed lines show 
^®Ni-powered cases. 


Fig. 4. — SN light curve of pulsar-driven SN models around the 
peak. The gray lines with arrows indicate the observed range of 
SNe Ibc and BL-SN Ibc. 


ejecta velocity is 14j ~ 20,000 km s“^, which is com¬ 
patible with the observed BL-SNe Ibc. An interesting 
possibility is that SL-SNe Ic and BL-SNe are connected 
sequences, and the main difference is the strength of the 
magnetic field. 

Although the pulsar-driven SN model can explain the 
peak light curve of BL-SNe Ibc, the radioactive decay of 
®®Ni has typically been considered as the main energy 
source, so as in the case of ordinary SNe Ibc. The peak 
luminosity powered by the ®®Ni decay can be roughly 
estimated as ~ x or 


r=®Ni 


4 X 10'^^ erg s ^ 


/ MseNi \ / Mej \ 

\0.1MqJ \5Mq) 


\10® cm ) \0.05 g“^ cm^ 


( 12 .) 


On the other hand, the observed bolometric luminosi¬ 
ties range from ^ j^q 42-43 gj.g g-i ^ 

10'^^ erg s“^ for BL-SN Ibc. The synthesized ®®Ni masses 
are estimated to be ^ 0. 05 — 0.8 Mp,, although the un- 
certeinties are large ('e.g.. lDrout et al.l2011l : ILvman et al.l 

[Ml . 

Fig. 13] shows several sample light curves. The blue 
dashed lines are the cases in which only ®®Ni decay is 
considered. The gray lines are the observed l ight curves 
of SNe Ibc and BL-SN Ibc (iDrout et al.ll201lll . Compar¬ 
ing Eqs. m and m, one sees that the pulsar-driven 
model may also mimic SN light curves, with the flux 
as dim as that of observed SN Ibc by considering a rela¬ 
tively large magnetic fields, Bdip ^ 5 x 10^® G. Note that 
a relatively slow rotation of Pi > 10 ms better explains 
ordinary SNe Ibc; if the spin is faster, the SN ejecta is 
inevitably accelerated up to a high velocity and the SN 
becomes brighter (see Eq. jH]). At this stage, one could 
speculate that some of the BL-SNe Ibc and SNe Ibc are 
also connected sequences. Both can be driven or aided 
by newborn pulsars with a magnetar-class dipole field 
and the difference is the spin. 



1 10 


Pj [ms] 

Fig. 5.— Contour plot showing properties of SN counter¬ 
part of fast-spinning strongly magnetized proto-NS formed with 
Mej = 2 Mq, Mssni = 0-05 Esn = 1 x 10®^ erg and 

K'y = 0.05 g~^ cm^. The color with solid lines shows the peak 
absolute magnitude, the dotted-dashed lines show the decline rate 
of the light curve, Mr^s = 0.3 and 1, and the dotted line shows 
the contour of Sk = 1 x 10®^ erg. The parameter region broadly 
consistent with the observed SN Ibc is indicated. 

3.1. Optical Constraints on Pi and B^ip 

In addition to the peak luminosities we discussed 
above, rising and decaying timescales of SN light curves 
can be used to constrain physical parameters of underly¬ 
ing proto-NSs. Fig. Hjfocuses on the raising and early de¬ 
cline of light curves. The gray li nes indicate the ob served 
range of SN Ibc and BL-SN Ibc (jDrout et al.ll20Illl . The 
decline rate is in the range of 0.3 < < 1, where 

Mr_i 5 is Mr, — A/R_niax at 15 days after the peak. The 
thick solid red line shows a pulsar-driven case broadly 
consistent with the observed SNe Ibc. The evolution of 
a light curve becomes wider when the poloidal field is 
smaller (dash line) because the energy injection rate de¬ 
clines more slowly. Also, a larger ejecta mass case (dot¬ 
ted dash line) gives a slow light curve because the photon 
diffusion time becomes longer. 

In Figs. [5| and [O] we show in what parameter range 


Peak R-Band AB Magnitude 
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Fig. 6. — Same as Fig. [5] but with Mgj = 5 Mq, Msejsji “ 
0.1 Mq, and £sn = 3 x 10^^ erg. Parameter regions broadly con¬ 
sistent with the observed SN Ibc, BL-SN Ibc, and SL-SN Ic are 
indicated. 


the pulsar-driven SN model can explain the observed op¬ 
tical emission from SE SNe. Fig. [Mel) corresponds to 
relatively low (high) ejecta mass, Mej = 2 Mq (5 Mq). 
In both cases, the ®®Ni mass and SN explosion energy 
is moderate and the SN emission is predominantly pow¬ 
ered by the magnetized wind except for the bottom right 
Conner of the panels. The boundary of the Ni domi¬ 
nated region is shown by the solid white line. SNe Ibc 
with Mu^max —(17 — 18) and ^0.3—1 can 

be explained by the pulsar-driven SN model in the top 
right Conner of the panels, Tj > 10 ms and > 

5 X 10^"^ G. Note that proto-NSs with relatively weak 
poloidal fields cannot hide in SNe Ibc since the light 
curves become slower than the observed ones. BL-SNe 
Ibc with Mu^niax —(18 — 19) and I4j ^ 20, 000 km s“^ 
can be explained by a larger-mass case, Mej > 5 Mq, 
with Pi about a few ms and i?dip ^ 5 x 10^“^ G, in which 
the kinetic energy is also mainly provided by the magne¬ 
tized wind. SL-SNe Ic also prefer relatively large ejecta 
mass cases since their light curves are relatively slow, 
the decrease in mag nitude 40 days after peak is < 1.5; 
(jQuimbv et al. 1170131 ). The best fitting parameter range is 
Pi less than a few ms and i?dip ^ 10^^ G. 

The possibility that a significant fraction of SE-SNe 
are driven by nascent pulsars is interesting in view of 
the connection among G RBs, SL-SNe and BL-SNe (see 
also [Metzger et al.l [2015[1 . It is also of interest in view 
of the connection to Galactic magnetars in the dynamo 
hypothesis. 


3.2. Late-time behavior 

As shown above, peak optical light curves of SNe Ibc 
and BL-SNe Ibc can be broadly explained by the pulsar- 
driven model with the appropriate choice of Pi and i?dip- 
On the other hand, these SNe have been considered to 
be mainly powered by ^®Ni decay. The parameter degen¬ 
eracy between Pi, Bp, and MseNi cannot be solved only 
from the peak optical light curves. One promising way is 
to use late-time spectroscopy. Indeed, in some cases, the 
®®Ni masses are independently determined by observing 
Fe line emissions in the Go decay phase (> 100 days), 
and are consistent with the values obtained from the 


peak optical light curves. However, such observations 
are challenging since the line emissions are typically very 
faint. Also, there are still significant uncertainties in the 
line transfer calculation. Another possible way to solve 
the parameter degeneracy is to use the late-time optical 
photometry from SNe, which can provide an independent 
constraint on the ®®Ni mass. However, it is known that 
the late-time light curves of SE SNe are heterogeneous 
and difficult to fit consistently with th e optical peak 
using a simple ^®Ni-decay model (e.g., iWheeler et al.l 
120151) . Also, pulsar-driven models could reproduce the 
light curves. Note that our simple model of calculating 
optical light curves becomes less reliable after the late 
decline phase or early nebular phase ~ 20 days after the 
peaks. More detailed theoretical calculations of late-time 
optical emission are necessary. 

4. MULTI-MESSENGER TESTS 

Because of additional parameters in the pulsar-driven 
SN model, optical light curves alone may not be used to 
distinguish the model from the other competing models. 
Multi-messenger approaches are useful to break param¬ 
eter degeneracies, to test the pulsar-driven scenario for 
SE SNe from ordinary SN Ibc to BL-SN Ibc and SLSN Ic 
and also the Galactic magnetar connection to SE SNe. A 
unique signature of ne wborn pulsar engines is the PWN 
emission in X rays (e.g., jPerna et al.ll2008l : iMetzger et al. | 

2014 ; [Murase et al.l[26l5[) a nd gamma ravs ([Kotera et ah 

2013 : [Murase et al.[ [2015[) . Although the dissipation 
mechanism of the magnetized wind is still controversial, 
a most likely outcome is an injection of ultra-relativistic 
electrons, which triggers leptonic pair cascades mediated 
via synchrotron emission and (inverse) Gompton scat¬ 
tering. The synthesized nebular emissions are entirely 
down-scattered into the thermal bath in the earlier phase 
of the ejecta expansion, but start to escape the ejecta at 
a later time. By observing such broadband nebular emis¬ 
sions in soft X-ray, hard X-ray, and gamma-ray bands, it 
is possible to put independent constraints on the physi¬ 
cal parameters of underlying NSs. Such signals can also 
probe the particle acceleration in embryonic PWNe. 

Moreover, fast-spinning strongly magnetized proto- 
NSs are possible sources of new messengers. In 
general, fast-spinning proto-NSs are unstable to non- 
axisymmetric perturbations and can evol ve into a plausi¬ 
ble configuration f or emitting GWs fe.g.. [Kokkotasl[2008l : 
[Bartos et al.l[2013l) . The GW frequency is / ~ 100 Hz — 
1 kHz, which coincides with the target frequency range of 
ground-based interferometers. In principle, the detection 
of such GWs can be used to determine physical param¬ 
eters of newborn pulsars, e.g., the rotation period and 
deformation rate. Neutrinos are also a powerful messen¬ 
gers. In addition to multi-MeV thermal neutrinos from 
proto-NSs, some hadron acceleration processes can occur 
in the magnetized wind or jet, and the ener gy dissipation 
resul t s in GeV to EeV neutr i no emissions (iMurase et al.l 
12014120091 : iFang et nil20l4 iLemoine et al.ll2015h . Such 
high-energy neutrinos can be a probe of the physics in 
strongly magnetized winds. 

4.1. High-energy X-ray and gamma-ray emission 

Non-thermal emission from PWNe can probe underly¬ 
ing newborn pulsar engines. The injection spectrum is 
a hard power law, dN^/dE-f oc E~^ with s ^ 1.5 — 2.5 
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Fig. 7.— Hard X-ray (30 — 80 keV) light curves of pulsar-driven 
SN model with typical parameters for SN Ibc, BL-SN Ibc, and SL- 
SN Ic labeled in Figs. |3and|6] The 3a detection threshold by 50 
ks observation using NuSTAR from 20 Mpc and 300 Mpc is shown. 


from soft X rays to GeV-TeV gamma rays (|Murase et al.l 
[Ml see also SecjC]). The light curve depends on the 
spin-down of the underlying NS. 

Here, we focus on the hard X-ray counterpart, where 
the Compton scattering is the main interaction process 
inside the SN ejecta and our theoretical calculation is 
most robust. We diyuss the detectability using NuS¬ 
TAR (|Harrison et al.l[2013f) . which operates in the band 
from 3 to 79 keV. Hard X rays can be also produced 
in the ®®Ni-powered model; the gamma-rays produced 
by the ^®Ni decay into ®®Co and ^®Fe with > MeV 
are successively Compton-scattered down to lower ener¬ 
gies. However, such hard X rays begin to be suppressed 
once the SN ejecta becomes Compton thin, while the 
PWN emission rises at the same moment. Moreover, the 
spectrum from Ni decay have a lower energy cutoff at 
^ 100 keV le.g.. lMaedall2006l l. and can be distinguish¬ 
able from the PWN spectrum. 

We should note that the PWN emission in other energy 
bands can be a useful counterpart too. If the energy in¬ 
jection by the pulsar wind is large enough and the ejecta 
mass is relatively small, the ionization break occ urs and 
soft-X-ray can be observed (jMetzger et al.ll2?iHl . Such 
a situation is promising in SLSN-Ic, but not guaranteed 
in BL-SN Ibc and SN Ibc. Also, the GeV gamma-ray 
counterpart can be better than X rays as a probe since 
it does not depend on the ionization state. Gamma-ray 
detection is possible for nearby o bjects up to < 10 Mpc 
as shown in iMurase et al.l (j2015[ l. 

Fig. [3 shows light curve examples of pulsar-driven SN 
model in the hard-X-ray band (30 — 80 keV). The thick, 
dash, and dotted-dashed lines correspond to the cases 
for which the SN counterpart is consistent with SN Ibc, 
BL-SN Ibc, and SL-SN Ic labeled in Figs. [5] and El re¬ 
spectively. We also show the 3 (t detection threshold us¬ 
ing NuSTAR with 50 ks observation. The hard-X-ray 
counterpart can be detectable for SN Ibc and BL-SN 
Ibc at < 30 — 50 Mpc and for SL-SN Ic at < 1 Gpc 
{z < 0.2). The anticipated detectable event rate is 
^ 1 yr“^ sky“^ for SN Ibc, BL-SN Ibc, and SLSN-Ic 
(see Fig. [1]). A follow-up observation needs to be under¬ 
taken ~ 50 — 100 days after the explosion. 


Mei = 2M,„„,E,„=lxlO^‘erg 



1 10 
Pj [ms] 


Fig. 8. — Contour plot showing properties of hard X-ray counter¬ 
part of fast-spinning strongly magnetized proto-NS formation with 
Mej = 2 Mq and £’sn = 1 x 10^^ erg. The color with solid lines 
corresponds to the 3 (j detection horizon using NuSTAR with 50 ks 
observations at around the peak time shown in the dotted-dashed 
lines. 

The raising time of the hard X-ray counterpart can be 
roughly estimated from the condition AlcompFcomp^ ~ 1, 
which gives 


tx ^ 80 days 


Mej 
5 Mo 


( v-j V7 

\10®cms“^/ \50keVy 

.03) 

In the above estimate, we approximate the inelasticity 
of Compton scattering as ATcomp ~ E-^/meC^{^ I). At 
t > tx, the PWN emission can be directly observed, i.e.. 


fe: 


1. Here, fesciE-y) is the fraction of the PWN 


emission escaping from the SN ejecta (see Eq. ICI7I) . The 
peak luminosity is roughly given by ~ Lem,i/7^b x 
where the bolometric factor becomes TZh 
10 — 20 around the hard X-ray raising time. From Eq. 
(l9|), the ratio between the SN and hard X-ray luminosity 
is given by 

Tx 1 Ki 


rPsr 

-^sn 


7^bi^c0mp C 


^0.007 


lb) 


Eei 


10® cm s“i / V 50 keV 


E^ 


- 1/2 

(14) 


These r esults are cons i stent with more detailed calcula¬ 
tions bv iMurase et al.l (j2015fl . 

The detectability of the hard X-ray counterpart is 
shown for more general cases in Figs. El and El for which 
we use the same parameter set as in Figs. [5] and El- The 
color contour with solid lines shows the detection hori¬ 
zon using NuSTAR, and the dotted-dashed lines show 
the emission raising time. The hard X-ray counterpart 
is a promising signature of identifying newborn pulsar 
engines. 


4.2. Gravitational wave emission 

In general, fast-spinning proto-NSs are unstable to 
non-axisymmetric perturbations and can evo lve into a 
plausible configura.tion fo r emitting GWs (e.g.. lKokkotasl 
120081 iBartos et aLll20I3ll . If the energy loss through the 
GWs is significant, then it might suppress the electro¬ 
magnetic counterparts. Thus, it is useful to consistently 


NuSTAR Detection Horizon [loglO Mpc] 
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Fig. 9.— Same as Fig. [S] but with Mej = 5 Mq and Ssn = 
3 X 10®^ erg. 

model GW spin-down together with the electromagnetic 
spin-down. We focus on t he GW emiss i on due to magnet¬ 
ically deformed rotati on (|Cutlei]l2002t iStella et al.ll2005t 
iDairOsso et al.l 120091) , which is an interesting channel 
especially in terms of magnetar formation. A proto-NS 
with a strong inner toroidal field Bt is deformed by the 
magnetic pinch effect. In general, the axis of the defor¬ 
mation is different from that of rotation and the proto- 
NS starts to precess. The tilt angle of the precession in¬ 
creases secularly due to the bulk viscosity and the proto- 
NS evolves into a non-axisymmetric rotating body, which 
is a plausible configuration for the GW emission (see the 
Appendix). The GW form is parameterized by Pi, Bp 
and Bt (or, the deformation rate, ec). 

First, let us argue the detectability of the GW coun¬ 
terpart under the competition with the electromag¬ 
netic spin-down. For a given spin-down timescale, i.e., 
P{t) and GW luminosity Lgw(t), the signal-to-noise ra¬ 
tio (S/N) of the expected GW averaged over all possi¬ 
ble orientations of source an d detector can be estimated 
as (|Owen fc Lindblomll2002ll : 



4G 

f dJ 

UJ 

IOttc®!!^ j 

fSM) 


Here, / = 2/P is the GW frequency, dJ = dLg^^/{2 'it/P), 
and S'h(/) is the one-sided power spectral density of de¬ 
tector noise. We note that the {S/N) calculated from Eq. 
(fT5|) is roughly equal to that obtained by the matched- 
filtering analysis. If the excess-power search is imple¬ 
mented, which is more appropriate for this type of GW, 
the anticipated (S/ N) would become smaller by a fac- 
tor of a few to ~ 10 (|Thrane et al.l[2Mll : iPiro fc Thrand 
nnn). We use an anticipated sensitivity curve of Ad¬ 
vanced LIGO which is for the optimal direction of the 
detector and the angle-averaged sensitivity is smaller by 
a factor of ~ 2/3. On the other hand, by combing other 
detectors, e.g.. Advanced Virgo and KAGRA, the sensi¬ 
tivity effectively increases at most by a factor of ^ -x/S. 

Fig. [TUI shows the detection horizon of the GW coun¬ 
terpart. Each panel shows a different deformation rate; 
Bt = 1.0 X IQi® G (ec = 3.3 x IQ-^; top), Bt = 

https://dcc.ligo.org/LIGO-T0900288/public 
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Fig. 10.— Detectability of the GW counterpart of magnetically 
deformed proto-NS. The solid-line contour with color shows the 
detection horizon of c/l = 5,10,15, and 20 Mpc with (S/N) = 8 
by using Advanced LIGO. The dotted-dash-line contour shows the 
spin-down timescale of the proto-NS, tgd = 0-1? 1^)0 days. 

2.0 X 10^® G (ec = 1.3 x 10“^; middle), and Bt = 
3.0 X 10^® G (ec = 3.3 x 10“®; bottom). Such a defor¬ 
mation rate is recently inferred for som e Galactic magne- 
tars f r om th e X-ray timing observation (IMakishima et all 
12014 120151) . The solid-line contour shows dh = 
5,10,15, and 20 Mpc with (S/N) = 8, which is the 
stan dard threshold value for compact binar y merg¬ 
ers (|LIGO Scientific Gollaboration et al.l l2013|) . The 
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Fig. 11.— SN light curves of pulsar-driven SN models with 
millisecond rotation and different toroidal magnetic field strength, 
St = (1 — 3) X 10^^ G. We fix other parameters as S^ip = 2 x 
10^^ G, Pi = 1 ms, Mej = 2 M 0 , and Kt = 0.05 g“^ cm^. 

dotted-dashed line contour shows the spin-down 
timescale of the proto-NS, tgd = (1/t^™ + = 

0.1,1,10, and 100 days. In Fig. [10] we shut off the 
GW spindwon for a larger toroidal field following Eq. 
(|A5p . The spin-down timescale via GW emission can be 
roughly given by 



From Eqs. OT and (HU, the GW spin-down dominates 
when 

.o>3xl0-< (17, 

In general, the (S/N) becomes larger for a smaller dipole 
field because the competitive electromagnetic spin-down 
becomes irrelevant and for a faster rotation because the 
intrinsic energy budget becomes larger. In principle, the 
GW can be detectable up to the Virgo cluster, du = 
16.5 Mpc for Pi less than a few ms, B^ip less than a few x 
1013 G, and Bt>2x lO^^ G (ec > lO"^). 

With a sufficiently large (S/N), physical parame¬ 
ters like Pq and Bt (or ec) can be determined from 
the GWs. The determinatio n acc uracies can be at 
most APi/Pi Ri {S/N)~^/^Ncyc and ABt/Bt ^ 
Aec/ec « iS/N)-\ Here, 

6 X 103(eG/10“3)-i(P;/ins)3/2 corresponds to the num¬ 
ber of GW cycles in the spin-down timescale. Thus, if 
this type of GW is detected, the rotation period of proto- 
NS could be determined with a sufficient accuracy. 

Next, let us discuss the effect of GW spin-down on 
the electromagnetic counterpart. Fig. [TT] shows sev¬ 
eral sample light curves of the pulsar-driven SN model 
with millisecond rotation and different toroidal magnetic 
field strength. We set Hdip = 2 x lO^^ G, P, = 1 ms, 
Mej = 2 Mq, and Kt = 0.05 g“i cm^. The SN emission 
becomes dimmer for a stronger GW spin-down, from the 
SL-SN class to the ordinary SN Ibc class. A broader pa¬ 
rameter region is investigated in Figs. 1121 and [T51 where 
we assume the same parameter set as in Figs. |5| and |6| 
except for Bt = 3.0 x lO^® G (eq = 3.0 x 10 “ 3 ). Gom- 


Mej = 2 M,„„, M^i = 0.05 M,„,„ E,„ = 1 X 10^‘ erg, Eq = 3.0 x lO'^ 



1 10 
Pj [ms] 


Fig. 12.— Contour plot showing properties of the SN coun¬ 
terpart of fast-spinning strongly magnetized proto-NS formation 
with Mej = 2 Mq, Mse^i = 0-05 Mq, fan = 1 X 10®^ erg, and 
£0 = 3.0 X 10“®. The color with solid lines shows the peak abso¬ 
lute magnitude, the dotted-dashed lines show the decline rate of 
the light curve, Mr^s = 0.3 and 1.0, and the dotted line shows 
the contour of fx = 1 x 10 ®^ erg. 



Fig. 13. — Same as Fig. [S] but with Mej = 5 Mq, Mssni = 
0.1 Mq, fan = 3 X 10®l erg. 

paring with Figs. [SjandlHl the peak magnitude becomes 
significantly smaller in the parameter region where the 
GW spin-down is relevant. As a result, for a relatively 
small ejecta mass case (Fig. |T2|), the bottom left Conner 
{Pi about a few ms and Pdip less than a few x 10^3 Qj 
becomes consistent with ordinary SN Ibc. As for a rel¬ 
atively large ejecta mass case (Fig. |T3|), the parame¬ 
ter region compatible to SL-SN Ic disappears due to the 
GW spin-down. For a larger poloidal magnetic fields, 
^dip ^ 10^“^ G, the effect of GW spin-down is not notice¬ 
able. 

Fig. [T4| shows the hard X-ray counterpart of the 
SN shown in Fig. |TT| The hard X-ray counterpart is 
also more suppressed for a larger GW spin-down. A 
broader parameter space is investigated in Fig. |T5] for 
Mej = 2 Mq cases. The hard X-ray counterpart of 
pulsar-driven SN Ibc with strong GW spin-down can be 
detectable even from ^ 100 Mpc by follow-up observa¬ 
tions ^ 50 — 100 days after the explosion using NuS- 
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Fig. 14.— Hard X-ray (30 — 80 keV) light curves of pulsar- 
driven SN models with millisecond rotation and different toroidal 
magnetic field strength, Bt = (1 — 3) X 10^® G. We fix other 
parameters as fJuip = 2 X 10^® G, Pi = 1 ms, and Mej = 2 Mq. 



Fig. 15. — Gontour plot showing properties of hard X-ray coun¬ 
terparts of fast-spinning strongly magnetized proto-NS formation 
with Mej = 2 M©, Mse^i = 0.05 M©, fan = 1 X 10®^ erg, and 
CQ = 3.0 X 10“®. The color with solid lines corresponds to the de¬ 
tection horizon using NuSTAR with 50 ks observations at around 
the peak time shown in the dotted-dashed lines. 

TAR. The maximum detectable event rate both by the 
second generation GW detector network and NuSTAR is 
^ 1 yr“^ sky“^. 


4.3. Neutrino emission 

We briefly overview possible neutrino counterparts 
from fast-spinning newborn NSs. Neutrinos can be used 
as a powerful probe of fast-rotating pulsars, because they 
can escape without significant attenuation. Multi-MeV 
thermal neutrinos fr om a nearby SN are detectable with 
Hyper-Kamiokande (I Abe et al.ll2011[l . In addition, if the 
proto-NS is rotating and magnetized, nucleons supplied 
by the neutrino-driven wind should be accelerated mag¬ 
netically, leading to quasi-thermal neutr ino emission in 
the GeV-TeV range (jMiirase et al.l[2?)T^ . Even higher- 
energy neutrinos can be p roduced depending on the mag- 


netic dissipation process ( 

Murase et alJl2009l: iFane et alJ 

l2014HLemoine et al.ll2015 

). If the pair-multiplicity is not 


as much large, PeV-EeV neutrinos, which can be pro¬ 


duced by pp and/or py interactions, can be detected up 
to ^ 10 Mpc. Note that high-energy neutrino emission 
is typically expected for SNe II when the S N shock be¬ 
comes collisionless after i ts shock breakout (iKatz et al.l 
120111 : iMiirase et al.ll2011l : FKashivama et al. I [201,3( 1. 

5. DISGUSSION 
5.1. Observational strategies 

Based on the results in the previous sections, we dis¬ 
cuss possible observational strategies for fast-spinning 
strongly magnetized newborn NSs in SE SNe. 

5.1.1. SNe Ibc 

The optical light curve of SN Ibc can be broadly consis¬ 
tent with the pulsar-driven SN model with Pj > 10 ms, 
.Bdip ^ 5 X 10^"^ G and Mej < 5 Mq. For such cases, 
NuSTAR can detect the hard X-ray PWN emission from 
dh ^ 30 Mpc. We encourage follow-up observations of 
^ 50 ks, ^ 50 — 100 days after the SN explosion. The 
anticipated SN Ibc rate within the NuSTAR detection 
horizon is about a few yr“^ (see Fig. [T]). 

If the inner toroidal magnetic field is as strong as 
Bt ^ 10^® G, the GW spin-down due to magnetically de¬ 
formed rotation can effectively suppress the pulsar wind, 
and the pulsar-driven SN emission can be consistent with 
SNe Ibc if Pi is about a few ms and Pdip 10^“^ G. Such 
GWs from SNe Ibc can be detectable up to < 20 Mpc 
by Advanced LIGO. NuSTAR can also detect the hard 
X-ray counterpart with a same follow-up observation as 
above. Note that the typical GW emission duration 
is ^ 0.1 days, which is much shorter than the rele¬ 
vant timescales of the SN and hard X-ray counterpart, 
> 10 days. In order to pin down the GW time win- 
dow, it will be useful to c ombine neutrino emission (e.g., 
iMurase et al.l I2009L1201411 or some electromagnetic pre¬ 
curso r signals (e.s.. iKistler et aO 120131: iNakauchi et al.l 

nnn). 

Multi-messenger detections or even non-detections 
from nearby SNe Ibc give important implications on 
the origin of Galactic magnetars. If the progenitors of 
Galactic magnetars are WRs and the formation rate is 
as high as Eq. (SD, ordinary SN Ibc is the promis¬ 
ing site of the magnetar formation. Detections of the 
hard X-ray counterpart can fix where and how Galac¬ 
tic magnetars are formed. On the other hand, non¬ 
detection may indicate a lower formation of magnetar 
rate than Eq. (|3]) or different formation pathway, e.g., the 
“fossil” scenario (lUsovlll99^lThompson fc Duncani[l993t 

iFerrario fc Wickramasinghell2006h . Moreover, the multi- 

messenger detections can strongly constrain the initial 
spin of proto-NS, which is a missing link in the massive 
star evolution; spins of ma ssive 0-type st ars and NS pul¬ 
sars have been measured (lFukudalll982D . but no obser¬ 
vational information about the in-between. Probing the 
proto-NS spin will shed light on the angular-momentum- 
transfer process during m assive star evolution , which is 
still highly uncertain ('e.g.. IMevnet et ani2011[l . 

5.1.2. BL-SNe Ibc 

The optical light curve of BL-SN Ibc can be broadly 
consistent with the pulsar-driven SN model with Pj less 
than a few ms, Pdip ^ 5 x 10^“^ G, and Mej > 5 Mq. 
For such cases, NuSTAR can detect the hard X-ray PWN 
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emission from cJl ^ 50 Mpc. It requires a follow-up ob¬ 
servation of ^ 50 ks, ^ 50 — 100 days after the SN explo¬ 
sion. The BL-SN Ibc rate within the NuSTAR detection 
horizon is ~ 1 yr“^ (see Fig. [1]). 

Detections of such hard X-ray counterparts will indi¬ 
cate that newborn pulsars play relevant roles in BL-SNe 
Ibc. Spin-down power of newborn pulsars can be as im¬ 
portant as ®®Ni decay for the SN emission. Also, the 
large kinetic energy inferred for BL-SN Ibc ejecta can be 
mainly provided by the spin-down power. The pulsar- 
driven SN model for BL-SN Ibc requires a magnetar-class 
field strength and Pi ^ ms. So, the hard X-ray observa¬ 
tion can probe the nascent stage of magnetars, namely 
indicating that the ms rotation is relevant for the mag¬ 
netic field amplification. We should note, however, that 
a Galactic magnetar r emnant is known to b e less ener¬ 
getic, £k ^ 10^^ erg (|Vink fc Kuioeii [2006H . Also, the 
observed BL-SN Ibc rate is smaller than the magnetar 
formation rate, thus BL-SNe Ibc could explain only a 
minor abundance of Galactic magnetars. 

5 . 1 . 3 . SL-SNe Ic 

The pulsar-driven SN model can explain the optical 
light curve of SL-SN Ic with Pi less than a few ms 
and i?dip ^ 10^^ G. For such cases, the hard X- 
ray counterpart can be detectable from ^ Gpc 
{z < 0.2) using NuSTAR with ^ 50 ks observations. 
A follow-up observation ^ 100 day after the SN ex¬ 
plosion is required. Gur rent optical tra nsient surveys 
like P TF (|Law et al.ll2009ll. Pan-STARRs (|HodaDD et al.l 
1200411 . and AS ASSN (iShaopee et al.l[20T^ can detect the 
optical counterpart with a rate larger than a few yr“^ 
from within the NuSTAR detection horizon. 

Detections of the hard X-ray counterpart can strongly 
support the pulsar-driven SN model for SL-SNe Ic. We 
should note that early PWN emissions could be also ob¬ 
served in soft X ravs (IMetzger et al.11201411. GeV-TeV 
gamma rays (|Kotera et al.l 120131 IMurase et al.l 1201511 . 
and radio, but the detectability will be more sensitive 
depends on, e.g., the ionization fraction of the ejecta and 
the lepton acceleration in the PWN. 

5 . 2 . Impacts of Simplifications 

Our semi-analytic model includes several simplifica¬ 
tions, which needs to be refined for more detailed com¬ 
parison with observation. The SN light curve is cal¬ 
culated based on the one-zone approximation. Includ¬ 
ing the multi-dimensional effects is crucial to obtain 
multi-band light curve more precisely. We effectively 
fix the opacity of the elastic scattering and photoelec¬ 
tric absorption separately, but these quantities need to 
be determined consistently in a time-dependent man¬ 
ner. More detailed radiation-transfer calculations tak¬ 
ing into account the ionization degree of metals are re¬ 
quired. Our treatment basically overestimates the pho¬ 
toelectric absorption, and thus the SN counterpart, and 
underestimates the soft X-ray counterpart, although the 
effect is minor for the SN counterpart. As for the in¬ 
jection spectrum from PWNe, we use a simple broken 
po wer law motivated b y detailed numerical calculations 
bv IMurase et al.l (I2015D . However, the shape of the spec¬ 


trum in general changes with time depending on the 
Compton parameter Y in the nebula. Our treatment 
overestimates the gamma-ray flux once Y becomes small. 

6. SUMMARY 

To test the pulsar-driven supernova (SN) scenario for 
stripped-envelope (SE) SNe from broad-line SNe Ibc 
to super-luminouse SNe Ic and ordinary SN Ibc, and 
also the Galactic magnetar connection to SE-SNe, we 
calculate multi-messenger counterpart of fast-spinning 
strongly magnetized proto-neutron star (NS) formation 
in massive collapses with SE. The SN emission powered 
by pulsar spin-down and ®®Ni-decay, early pulsar-wind- 
nebular (PWN) emission, and gravitational wave (GW) 
spin-down are consistently modeled. 

We show that the peak light curves of all types of ob¬ 
served SE-SNe, can be broadly explained by the pulsar- 
driven SN model; Pi > 10 ms and B^ip ^ 5 x 10^^ G for 
SN Ibc, Pi less than a few ms and Bup ^ 5 x 10^"'^ G for 
BL-SN Ibc, and Pi less than a few ms and Hdip > 10^^ G 
for SL-SN Ic. The latter two cases prefer more massive 
progenitors. 

For all cases, the early PWN emission especially in the 
hard X-ray band can be the smoking-gun signal of an 
underlying newborn pulsar engine, detectable by follow¬ 
up observations using NuSTAR ^ 50—100 days after the 
explosion. The hard X-ray detection horizon is < 30 Mpc 
for SN Ibc and BL-SN Ibc, and < 1 Gpc for SL-SN Ic, 
and the potential detection rates are ~ 1 yr“^ sky“^. 

If the inner toroidal magnetic field is as strong as 
Bt 10^^ G, the GW spin-down due to magnetically de¬ 
formed rotation can be relevant especially for the cases 
with Pi about a few ms and Pdip ^ 10^^ G. The GW 
counterpart can be detectable up to < 20 Mpc by the 
second generation GW detector network. When the GW 
spin-down is strong enough to be detected, the pulsar- 
driven SN cannot be as bright as SL-SN Ic; instead, mil¬ 
lisecond proto-NSs with Pdip a few x 10^^ G result in 
ordinary SNe Ibc. 

When this work was completed, we became aware of 
the related, independent work, arXiv:1508.02712, While 
a part of the interest is shared, our work has more 
focus on the scenarios explaining both broad-lined su¬ 
pernovae and superluminous supernovae, and multi¬ 
messenger prospects. 
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APPENDIX 


In this Appendix, we show the semi-analytic model we use for calculating the multi-messenger counterpart of fast¬ 
spinning strongly magnetized NS formation in massive star core collapse. 

SPIN-DOWN 

The spin-down of proto-NS is calculated from (jOstriker fc Gunnlll969H 

- ^ = Lem + Lgw. (Al) 

The first term in the right-hand side is the magnetic spin-down power; 

iem = ^^^^^J^(l + Csin2x^), (A2) 

where /i = /2 is the magnetic moment, P is the rotation period, x^^ is the angle between the magnetic and 

rotat i on axis, and C* ^ 1 is a pr e-factor. Eq. (IA2I) is o btained by force-free simulations (|Gruzinovl 120051 : iSoitkovsk^ 
120061: iTchekhovskov et al.l 1201311 . Note that Eq. (IA2I) is larger than the classical dipole spin-down luminosity by 
3(1 -f Csin^ x^)/2 sin^ > 1. We assume that the magnetized wind is isotropic for simplicity. These approxima¬ 
tions are not valid within the Kelvin-Helmholtz timescale ^ 100 s wh ere the baryon loading on the magnetized 
wind via neutrino-driven wind from the proto-NS surface is relevant (e.g., iThompson et al.ll2004f) . Also, for an ex¬ 
tremely strong poloidal field, i?dip > 10^® G, ms-proto-NSs spin down within the KH timescale. In such cases, the 
magnetized wind could punch out the pro genitor star as a bi-polar jet collimated by the anisotropic stress and the 
hoop stress (iBucciantini et al.ll2007L 1200811 . We here only consider a longer timescale t ^ tKH.i/ and a poloidal field 
Ldip < 1015 G. 

The second term in the right-hand side of Eq. (EU represents the GW energy loss; 

2GieGl)H2n/Pr . 2 2 ^ 

Lgw = --g-sm Xe(,(l-h 15sm Xej,), (A3) 

0 C'^ 


where ec = AI/I is the deformation r ate, P' is the pattern period, and Xeq angle between the deformation 


and rotation axis (IGutler fc Jones 


iStella et 1111120051: iDaH’Osso et al 


2001). In this paper, we consider the magnetically deformed rotation (IGutledl2002l : 


2009( 1. which is an interesting channel especially in terms of magnetar formation. 


Once inner tor oidal magnet i c fields amplified up to a magnetar value, the proto-NS becomes oblate by the magnetic 
pinch (see, e.g.. lGutieill2002l : iKiuchi fc Yoshidall2008HGualtieri et al.ll201lll . The deformation rate can be estimated as 


15 £b 
4 |W| 


2.5 X 10 


-A- 


Bt 


\ 1016 G 


(A4) 


Here, |IF| ~ M^sC^ x 0.6C/(1 — 0.5C) ^ 4.4 x lO^^ erg is the gravitational energy of the proto-NS with compactness 
parameter, C = GM^s/RnsC^ ^ 0.17 (iLattimer fc Prakashll200lll . In gener al, the deformation axis is not coincide 
with the rotation axis, and the proto-NS processes around the rota tion axis (IMe stel fc Takhai([T972ll . The tilt angle 
of the precession can increase secularly due to the bulk viscosity (|DairOsso et al.l 12009( 1. and the proto-NS evolves 
into a prolate shape, which is a plausible configura tion for the GW emission {xp. = Xea — P = P')- Recently, 
a processing motion driven by deformatio n as Eq. (1A4D was inferred for a galactic magnetar from the X-ray timing 
observation (|Makishima et all (2014 (2015(1 . This GW emission can only occur if the viscous dumping timescale of 
the NS precession is sho rter than the competitive magnetic braking timescale. This condition can be described 
as (|DairOsso et al.[[2009(l 


Bt < 2.4 X IO 15 G 





Rdip 

1014 G 


-2 

-f 1 


1/2 


(A5) 


DYNAMICS 


We here describe a simple model for dynamics of the SN ejecta and the resulting 
radius of the ejecta evolves as 


dt 


= K 


ej- 


electromagnetic emission. 


The 

(Bl) 


We assume the density structure of the SN ejecta as 


3 - (5 Mej 



(B2) 


We take 5 = 1 as a fiducial value for the index, and thus a dominate fraction of mass resides at around R ~ i?ej • Without 
significant energy injection after the explosion, the expansion velocity is almost constant, ~ 14j,i = (2Lsn/Mej)i/^ ~ 
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10,000 km s-i (Ssn/lO®^ erg)i/ 2 (Mej/M 0 )-i /2 cm s ^ up to the Sedov radius. On the other hand, when a fast-spinning 
proto-NS exists, the ejecta is accelerated by the magnetized wind. 


d£K 

dt 


tdyn 


(B3) 


where fx ~ Mejt4j/2 is the kinetic energy, Sint is the total internal energy, and 

Rei 

^dyn = 


(B4) 


is the dynamical timescale of the ejecta. The energy injection from the underlying pulsar occurs at the shock between 
pulsar wind and the SN ejecta. The radius of the shocked wind region evolves as 


dR^^ 

dt 




(B5) 


Here Vnh is the expansion velocity, which is obtained from the pressure equilibrium, J Lem x min[l, 
6PejV;b/7, or _ 


Hnb 


7 / Lem X min[l,T.^''t4j/c]dt 


6(3-(5) 


Me 


f 

\Rw 


3-S 


T^^Vej/c]dt/ATTRl 


(B6) 


The factor min[l, n^-'Mj/c] represents the fraction of the spin-down luminosity deposited in the SN ejecta (see Eq. IC4I 
for the definition of t^^). If i?w > Rej, we set i?w ~ -^ej- 


ELECTROMAGNETIC EMISSION 


The time evolution of the internal energy in the SN ejecta is described as 


dS'mt j 

~dr ^ ~ '■ 


- -1" /dep.emLem + /dep,®BNi-^®BNi + fdep,^^CoL^^Co- 

^dyn 


(Cl) 


The first and second terms on the right-hand side correspond to the energy loss via quasi-thermal SN emission and 
adiabatic expansion, whereas the third, foruth, and fifth terms correspond to the energy injection via magnetar wind, 
®®Ni and Co decay, respectively. 

The bolometric SN luminosity can be given by 



j ^ ^int 

j.ej ’ 

(C2) 

where 

+ej _ ej 4?ej 
‘'esc ^ T ^ ’ 

(C3) 

is the thermal photon escape time from the ejecta. 



Trp 

_ (3 - 6)KtM,; 

47ri?ej^ 

(C4) 


is the optical depth of the ejecta, /^erriu is the Thomson opacity, is the mean molecular weight per 

electron, rriu = 1-66 x 10“^"^ g is the atomic mass unit, and 0 < ^ < 1 is the effective ionization fraction. Since we are 
mainly interested in SE SNe, we take /ie = 2. In general, ^ depends on the temperature and composition and evolves 
with time. Here, for simplicity, we use fixed values ^ = 0.25 — 1, i.e., LTt 0.05 — 0.2 g“^ cm^, which is reasonable at 
around the optical peak of SE SNe. The temperature of the emission can be estimated as 

/ e. \ 

= (d|) 

where Vej ~ 47ri?gj/3 and a is the radiation constant. Note that the above method of calculating the SN emission is 
equivalent to the classical Arnett model (IArnettiri982[) with uniform ejecta temperature (|ChatzoDOulos et al.ll201^ . 

At the interface of the magnetized wind and SN ejecta, highly relativistic electrons are injected, which are further 
accelerated, e.g., at the shock or in the magnetic turbulence, and then rapidly cool via synchrotron emission and inverse 
Compton scattering. The scattered photons have very high energies so that they can produce pairs by two-photon 
annihi lation. The synthesized electron/positron is still energetic and produces another pair successively. (Murase et al.l 
(|2015f) numerically calculated the above electromagnetic cascade process by assuming the electron injection spectrum 
as 

die l(7e/7b)"® (7b < 7e < 7 m), 
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with <71 ^ 1 — 1.5, <72 ~ 2.5—3, 7 b 10^'® which is motivated by the observation of young PWNe (|Tanaka fc Takaharal 
I2010f) . The electron maximum energy can be estimated by equating the acceleration timescale face = ? 77 eWeC^/ei?c 
and synchrotron cooling timescale tsyn = 3mec/4crTt^B7e, he., 7 m ~ {QnejrjaTB)^^^. Here, the magnetic field energy 
density is given by 


Ub 


_ 3 / Lemdt 

Stt 47ri?^ 


(C6) 


with €b ^ 10 iMurase et al.l (|2015li found that the resultant nebula spectrum can be well approximated by 

broken power law; 




dN^ 

dE~f 


ee^em f {E,/ 

B-hE'^y^ \ {Ey/E^^y (^syn ^ ddy < e-y_niax)j 


a 


where TZh 2/(2 — qi) -|- ln(£:-,,_i„ax/E^yn) and Cg = 1 — es ~ 1- The low-energy part is dominated by the fast-cooling 
synchrotron emission from injection electrons with 7 e < 7b- The corresponding break photon energy is 


b 3 2^-8 


(C7) 


On the other hand, the high-energy part is mainly produced by inverse Compton scattering and successive pair cascade. 
The maximum energy is determined by the two-photon annihilation with the SN photons. 


^ 2kBEn ■ 


(C8) 


Injected non-thermal photons from the wind nebula are down-scattered, or absorbed during propagating through 
the SN ejecta. The main interaction channel depends on the photon energy; Bethe-Heitler (BH) pair production for 
E^ > 10 MeV, Compton scattering for 10 keV ^ Ej < lOMeV, photoelectric (bound-free) absorption for 10 eV < 
Bj ^ 10 keV, bound-bound and free-free absorption for lower energy bands. We calculate the energy deposition 
fraction of a photon as 

/dep = max[l, /d ep,sc “1“ /dep, ab]. (C9) 

The contribution from the Compton scattering is estimated as 

/dep.sc = 1 - (1 - , (CIO) 

where K^ rpn is the inelasticity of Compton scattering, (Tcomp is the Klein-Nishina cross section, and Tgomp is the optical 
depth fsee lMurase et al]l2015f) . On the other hand, the energy deposition fraction by the absorption processes can be 
expressed as 

/dep,ab = 1 6Xp( Tab), (^H) 

where 

Tab — Tgu -p Tpe, (C12) 

and tbh and Tpe are the optical depth of BH pair production and photoelectric absorption. We treat the BH pair 
production as an a bsorption process since the inelasticity is Kbb = 1 — 2l{E^lmf,(?) ~ 1 for E^/meC? ^ 2 (see 
IMurase et al.ll2015ll . The optical depth of the photoelectric absorption can be estimated as 


(3 - <5)ifpeMej 

where we use an approximate form of the opacity for oxygen-dominated ejecta, 

-3 


ETpe « 5.0 X C 


E^ 


10 keV 


? ^ cm^, 


(C13) 


(C14) 


Here 0 < C < 1 is a time-dependent factor determined from the effective ionization fraction. In this paper, we fix 
C = 0.5 for simplicity. Note that the fraction of the energy in the soft X-ray band (and lower energy bands) is always 
subdominant in our case. For example, the SN light curve around the peak does not change significantly depending 
on the details of the photoelectric absorption. The photoelectric absorption can be relevant for the SN light curve in 
the late phase (> 100 days after the explosion) and the soft X-ray light curve. 

The total energy deposition fraction of the magnetized wind is calculated by 


/d 


ep,em 


/ fdepiEry)Ey dEj 

jE/^dE, 


(C15) 
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Fig. 16.— Bolometric SN light curves obtained by this work (solid red line) and the Arnett model (green-dashed line). Both lines are 
calculated for the same input physical parameters; Sk = 1 x 10®^ erg, Mgj = 2 M@, = 0.1 Mq, and Kt = 0.1 g“^ cm^. 

where dN/dEj is the wind nebula spectrum and fdepiEj) is the energy deposition fraction of a photon with an energy 
Ej. On the other hand, the observed non-thermal nebula spectrum can be calculated as 


, dN -yobs 

dE^ 


/esc ^ E^ 


dN^ 

dE^' 


(C16) 


where 

/esc = {exp(-rx) + [1 - exp(-Tx)](l - (C17) 

X=comp,ab 

is the fraction of injected PWN emission directly escape from the SN ejecta. 

The energy injection rate from the ®®Ni decay is calculated as 

Tssxi = ikTsaxiesexi exp (-I , (C18) 

V tsejsliJ 

Lseco = M56Ni(e56Co - e56Ni)exp (-) , (C19) 

V “6Co / 

where Msexi is the ®®Ni mass, esexi = 3.9 x 10^° erg s“^ g“^, eco = 6.8 x 10® erg s“^ g“^, tsexi = 8.8 days, and 
^56 Co = 111-3 days. It is known that the energy deposition fraction from ®®Ni decay can be well approximated by 
/dep 56Ni = 1 ~ exp(—Teff) with TgS = (3 — /47ri?ej® and Kefi Ri 0.03 cm® g“^. In this paper, we instead use 

Eq. (|C9p to estimate the deposition fraction consistently with the wind dissipation, 


/dep,^^Ni(Co) 


y^/dep(g56xi(Co),»)'^S6Ni(Co),i£56xi(Co),l 

i 


y^7^56Ni(Co),ig56xi(Co),l 

i 


(C20) 


where £56Ni(Co).z and P56Ni(Co).z are the mean energy of the decay product and the decay probability. We consider 6 
decay channels for ^®Ni and 11 channels for ^®Co listed in iNadvozhinI (I1994D . We assume that all the energy of positron 
emission goes to the thermal bath. As we show in Fig. [161 our model can br oadl y reproduce SN light curves obtained 
by using the conventional Arnett model, which confirms the validity of Eq. (IC91) . 
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